DISCRETE MULTISCALE ANALYSIS: 
A BIATOMIC LATTICE SYSTEM * 



G.A. CASSATELLA CONTRA, + AND D. LEVI* 

Abstract. 

We discuss a discrete approach to the multiscale reductive perturbative method and apply it to a biatomic chain 
with a nonlinear interaction between the atoms. This system is important to describe the time evolution of localized 
solitonic excitations. 

We require that also the reduced equation be discrete. To do so coherently we need to discretize the time 
variable to be able to get asymptotic discrete waves and carry out a discrete multiscale expansion around them. 
Our resulting nonlinear equation will be a kind of discrete Nonlinear Schrodinger equation. If we make its continuum 
limit, we obtain the standard Nonlinear Schrodinger differential equation. 



1. INTRODUCTION. Nonlinear systems, and in particular nonlinear discrete systems, are 



gaining an increasing impact in modern science 33 



In 1955 Fermi, Pasta and Ulam (FPU) [Hj considered a unidimensional chain of atoms with 
nonlinear nearest neighbouring interaction to verify if nonlinearity could produce energy equipar- 
tition. Instead, they found recurrence, i.e. the motion of the chain for small energies was almost 



periodic [43]. To explain this result Kruskal and Zabusky found in 1965 42 a connection be- 
tween the FPU system and the Korteweg-De Vries equation (KdV), an equation introduced in 
fluid dynamics to describe one dimensional surface waves in the shallow water context [20] . By 
introducing the Inverse Scattering Transform, they where able to solve the Cauchy problem for 
the KdV equation 15 and to prove the existence of soliton solutions. 

In 1967 Toda 37 considered a dynamical system with exponential interaction, U(r) = er r + 
r — 1, the "Toda potential", whose small amplitude approximation gives the FPU system, and 
shares many of the integrability properties of the KdV equation. So the FPU system turns out to 
be an approximation of a discrete soliton model. 

Later more complicate atomic chains have been considered, as, for example, the biatomic one 
6,9,11,12, 16|[26] . These systems have various applications in physics and biology as, for example, 
in the study of ferroelectric perovskites, materials that, in certain crystallographic directions have 
an almost unidimensional frame, and in organic molecular chains. A biatomic chain of neighboring 
atoms A\ and A% is described by the discrete independent variable n and a continuous time t. 
However, the simplest nonlinear coupled lattice dynamical equations one can construct for this 
system are not solvable. Only special exact solutions may be found. 

Multiscale expansions 7,8,19 21,35,36 have proved to be important tools to find approximate 
solutions for many physical problems by reducing a given nonlinear partial differential equation 
to a simpler equation, which is often integrable 5 . Recently, few attempts to carry over this 
approach to partial difference equations have been proposed 2,10,22,23,32 . Almost all approaches 
considered contain some approximation, either based on physical or on mathematical reasoning as 
scaling transformations of the lattice provide a nonlocal result. In the following we prefer to stick 
to mathematical approximations as in this case it will be more evident what to do to improve the 
final result [l7] . 

In ref [9j a biatomic chain obtained as a first nonlinear approximation of a complex Lenard- 
Jones interaction between atoms has been considered. There the multiscale expansion of the 
continuous limit of the lattice model showed that the modulation of periodic solutions is governed 
by the Nonlinear Schrodinger differential Equation (NLSE). Here we consider the same model but 
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we are interested in carrying out the multiscale expansion on the lattice, i.e. we are looking for a 
lattice equation which in the asymptotic regime approximate the biatomic nonlinear lattice. To 
do so we need to discretize time to be able to allow for discrete asymptotic waves. If we keep a 
continuous time variable an asymptotic wave travelling on the lattice by necessity will be described 
by a continuous variable. So by necessity we go over to a differential system. 

Discretization of variables, besides representing an interesting problem in mathematical physics 
for its computerizability, it is also useful in itself. Measurements, for example, are based on 
sampling of physical variables such as space and time. It follows that physical models in which 
variables are defined on the lattice are easier to be compared with the real world we see in our 
measurements . 

In this work, we propose to continue the previous researches of biatomic chains considering 
both t and n as discrete variables. In particular, we shall assume, as these authors, that the system 
has a unharmonic cubic potential as in nature, potentials usually are non-symmetric. We shall 
thus apply a discrete multiscale reductive perturbative method to the model introduced by Campa 
et. al. [9] consisting of a biatomic chain with a nonlinear nearest neighbour interaction. 

In Section 2 we describe in detail the biatomic chain and write down the dynamical equations. 
Then in Section 3 we introduce some notions of discrete calculus and multiple scales defined on 
the lattice which we apply in Section 4 to the biatomic chain introduced in Section 2. In Section 
5 we analyze the resulting nonlinear discrete equation obtained and carry out its continuum limit. 
Finally, in Section 6 we draw some final conclusions. 

2. THE MODEL. We want to describe here a chain suitable to represent, for example, an 
a-helix channel, see Scott (1999) [33]. Our model consists of a biatomic chain formed by a sequence 
of pairs of neighboring atoms A\ and A 2 , with masses Mi and M 2 , respectively. Each pair, made 
of an atom of mass Mi and the following one of mass M 2 , can be considered as a " molec ule" . We 



denote by the index n the n molecule formed by the atom A\ and A2 (see Fig. 2.1 1. Let us 



A-i ^2 ^1 ^2 ^1 ^2 ^1 ^2 ^1 ^2 



y n 



Fig. 2.1. Pattern of a biatomic molecular chain in one dimension. The chain is formed by a sequence 
of pairs of neighboring atoms Ai and A2 . The displacements of the atoms of the molecule n are indicated 
with x n and y n . 

indicate with x n (t) and y n (t) the displacements of the atoms A\ and A2 belonging to the same 
molecule n. For each atom, we assume only nearest neighbourg interactions. Then, the total 
potential of the chain is given by 

U = ^2{Ux(y n - x n ) + U 2 (x n+ i - y n )}, (2.1) 

n 

where U\ is the intramolecular potential, between atoms belonging to the same molecule, and U2 
is the potential between different molecules. 
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Given a natural |3[|6[ |38| asymmetric potential with an absolute minimum in the equilibrium 
position as, for example, a Lenard-Jones potential, by taking the first terms of its Taylor expansion 
around the equilibrium position we can write the potentials U± and U 2 as 

Ui(r) = \k ir 2 + e -p lT \ U 2 {r) = l -k 2 r 2 + € -fi 2 r 3 , 

where k\ and k 2 are the harmonic constants, /3± and j3 2 are the cubic interaction constants and e 
is a small parameter which will play the role of the perturbative parameter. We assume that the 
interaction between atoms of the same site is stronger than that of atoms of different sites; thus 
k\>k 2 and |/?i|>|/32|- So, the Hamiltonian of our molecular chain turns out to be 

H = En{|[ M i^ + M 2V 2 n ] + \[ki{y n ~ x n ) 2 + k 2 (x n+1 - y n ) 2 ] 



2 

- 5 [A(2M - x n) 3 + P 2 {x n +i - y n ) 3 ]}, 



where x(t)= dx J-p and the equations of motion are 
3U 

M t x n = (2.2) 

= h(yn ~ x n ) - k 2 {x n - y n -i) + ePi(y n - x n ) 2 - e/3 2 (x n - y n -i) 2 , 
dU 

M 2 y n = — (2.3) 
9y n 

= -kl(y n - X n ) + k 2 (x n +l - Vn) ~ ePl(Vn - X n ) 2 + e/3 2 (x n +l - ynf- 



Eqs. (2.2 2.3) are a natural extension of the FPU model [14] to a biatomic system. 

3. MULTIPLE SCALES ON A LATTICE. Here we introduce the concepts necessary 
to extend the multiscale reductive perturbative approach introduced by Poincare [5] for the study 
of the asymptotic expansion of ordinary differential equations and extended by Taniuti to the 
reduction of partial differential equations 35 36 to the case of difference equations [17 24 32 . 

3.1. Lattices and functions defined on them. Given a lattice, we will denote by n the 
running index of the points separated by a constant spacing h. Thus to the lattice index n, we 
can associate a continuous variable x~nh defining the position of the points with respect to the 
origin, for convenience chosen to be with no loss of generality Xq = 0. 

If we introduce a small parameter e—N^ 1 , where N is a large integer positive number, we can 
define on the same lattice the slowly varying discrete variables rij (j=l,2,3,...) of constant spacing 
Hj and the continuous variables Xj (see Fig. 3.1) where 



n 



= N j nj, Xj = e j x. (3.1) 



If rij varies by one, n varies by N 3 , a number much larger than unity. For this reason, rij is 
a "slow variable" and provide an asymptotic behavior of the system. For each j there is a slow 
lattice variable corresponding to the slow index rij. rij will be an integer only if n is a common 
multiple of N ] ~ . 

Let us consider F n , a function of the discrete index n. An equation on the lattice is a functional 
relation which involves the function F at various lattice points, {F n +e}- In the case of the model 



considered before (2.2 2.3), I = ±1. We are interested to transform the system, defined on a 
lattice n, to the slowly varying lattices rij, providing the scales of the asymptotic behavior of the 
original system. This is equivalent to say that we are interested in transforming the system defined 
on x to the one with the slowly varying variables Xj. We can consider the function F n written in 
terms of the slowly varying lattice variables {rij}, with, for example, j = 1, 2, F n = f ni ,n 2 i an d we 
can carry out the e expansion of the function F n +e- 

Let us consider at the beginning the case of one slowly varying lattice n\, i.e. F n = f ni . As 
the shift operator T n acting on F n gives T n F n = F n +i, we have F n +e = T^F n - In order to extract 
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AO SO H| 120 100 nj(X|) 



16 H 2 24 32 n 2 [x2) 



Fig. 3.1. Rescaled lattices. 

the behaviour of the function F n+1 = F(x + h) on the new scales, let us carry out the Taylor 
expansion of F n+ i in powers of h. In such a case the shift operator can be expressed as an infinite 
order differential operator with respect to x, i.e. 

T n = eMhd x )=f^ (h ^. (3.2) 

fc=0 

Moreover, if we define a A operator as A« = (T n — l)/h, we have 

0, - Ml±^l, (3.3) 



and eq. (3.2) could be written as 



r ..fM«^. (3.4) 



fe=0 



Formulas (|3.3 3.4h are written in terms of A^ + \ However on the lattice we can define an infinite 



number of different difference operators which in the continuum limit, when h goes to zero, go 
over to the first order derivative. Among them it is important the asymmetric shift operator 
= jh{T n — T^ 1 ). In this case we have 

arcsinh(ft,Al s) ) -A (arcsinh(ft,Ai s) ))' £ 
d x = - h , -> T„ = ^ . (3.5) 

fc=0 



Introducing the slowly varying variable x\ and the corresponding lattice n\ in eq. (3.2), as d x 
ed Xl , we have 



Ti = e™* = = 1* = £ . (3.6) 

k=0 
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If we introduce more lattice variables, for example {rij}, with j = 1,2, then T n becomes 

T i = T ee T ee 2 = ^ {^hd Xl ) k g (£e 2 hd X2 y ^ ^ 



fc=0 



3=0 



Once we expand the operator d Xj in terms of shift operators we get an expression for F(n ± £) in 
terms of variations of f{n\,n2) with coefficients depending on e and I. 

As delta operators are linear combinations of shift operators, from eq. (3.10) it can be proved 
that for A = we have the following formula 



18 



24 



k\ 



(A^) k f ni =^-P(i,k)(AMyF n , 



(3.8) 



i — k 



where (A^) fc / ni is the fcth-difference of f ni respect to m, and the coefficients P(i,k) are given 
by P(i,j) = ^2 % a —jW a S°Gi t , where oj is the ratio of the increment in the lattice of variable n 
with respect to that of variable n\. In this case, taking into account equation (3.1), uj=N. The 
coefficients S" and G^ are the Stirling coefficients of the first kind and second kind, respectively. 
The result (3.8) can be inverted, providing: 



M 



(A^) k F n = J2~Q(i,k)(A^yf ni , 



(3.9) 



i—k 



where Q(i,j) is the same as P(i,j), but with w = N^ 1 = e. 

A general way to get these formulas is provided by the finite operator calculus [13,29,30 . The 
finite operator calculus prescribes the following formula [25] 



^ = £ 



k=0 



(£) k PkU) 

k\ 



(A„J fc , 



(3.10) 



where the functions Pk(j) are the unique basic sequence associated to the operator A ni , i.e. such 
that they satisfy the following conditions 



Po(ni) = 1, Pfc(0) = for all k > 0, 
A„ 1 p fc (n 1 ) = fcp fc _x(m). 

The basic sequences can be directly obtained by the transfer formulae: 

Pk(ni) = ) ni 

When A ni = A^ or A ni = An} , the basic sequences are: 



(3.11) 



pW(n 1 ) = ^n 1 ( 



Kd x , ' 



fe-i 



(xi)k = Xi(xi — h) ■ ■ ■ {xi — kh + h), 



(3.12) 



(3.13) 



p < C\n 1 ) = h k n 1 [ 



2hd x 



' k n 1 k - 1 = {2h) k G k { Xl ;-h,2h), 



where Gk(y,a,b) are the Gould polynomials 29 given by 



G k (y;a,b) 



y — ka 

y 

(y-ka)(b) k 



b > k 

(y — ka)(y — ka — b) ■ ■ ■ (y — ka — (k — 1)6). 



(3.14) 
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Let us also mention that for each A ni operator we can write from eq. (3.101 



1 ^ 1 r d j 



c^) 3 = ^2J^ UC^)*. 

fc=0 ' y 



(3.15) 



i.e. we can express the partial derivative as an infinite sum of differences whose coef ficien t s dep ends 
from the type of difference we are expanding into. In terms of A^ + \ from eqs. (3.10 3.131, eq. 



(3.6) reads 



oo 
fc=0 



while, in the symmetric difference case, it reads 



± n r n 



fe=0 



(2/ie) 
"fcT 



(3.16) 



(3.17) 



From eqs. ( 3. 16|3. 17 ) we get that any finite shift in the original equation will give rise to an 
expression in the slowly varying variables which involves an infinity of lattice points or, equivalently, 
contains differences at all orders of the function f ni . So to get a reduced equation on a finite number 
of points we need to cut the series by requiring that the function f ni be of finite order of variation. 
Let us introduce the following definition: 

Definition : The function f n is a slow varying function of order p if 



ap+7„ = 0. 



(3.18) 



Than we can prove the following Theorem: 

Theorem The function F n is a slow varying function of order p iff f ni is a slowly varying 
function of order p in its own variable, i.e. if A^j" f ni = 0. 

Proof. The proof of this theorem will be given in the case of A = A + , but it is easy to see 
that is valid for any delta operator. It is divided in two parts: 



a) Be f ni a slowly varying function of order p. From formula (3.9) it follows that 



AP +1 F n = 



OO 

£ 

i=p+l 



(P + 1)! 



(3.19) 



i.e. F n is also a slow function of order p. 



b) Be F n a slowly varying function of order p. From formula (3.8) it follows that 



A P+1 f 
a ni Jni 



oo 

£ 

t=p+i 



(P + 1)! 



P(i,p+l)A i n F n = 0, 



(3.20) 



i.e. f ni is also a slow function of order p. 



The expansion (3.17) can be performed in two steps: at first we write the shift operator in the 



n variable in terms of the derivatives with respect to %\ by formula (3.6) and then we expand the 



derivatives with respect to x\ in term of delta operators by formula (3.15). In doing so we will 



have formulas in derivatives which are valid for any delta operator. Moreover the first expansion 
has e dependent coefficients while the second will provide a finite number of terms only if we use 
the slow varying condition for the functions f ni ,n 2 - 

Let us now explicitate the first terms of eq. (3.17) for future use, at first in term of the 
derivatives and then in delta operators assuming that the function f ni ,n 2 is a slow function at 
most of order 2. At first we shall consider the case in which we have only one slow lattice, just 
the variable n\ is present and then we extend the result to the case of two slow lattices, n\ and ni 
and to partial lattices n and m. 
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3.1.1. F n =f nx =f( Xl ). FromEq. we get 



2^2 



h 2 e 



F n±1 = f( Xl ) ± hed x J( Xl ) + ^d 2 x J(xx) + 0(e 3 ). 



(3-21) 



(3.22) 



As from eq.(3.15l for p=2, d Xl =A ni and d Xl ={A ni ) , then equation (3.21) reads 

F n ±\ = f ni ± ^(/m+i _ fm-i) + 7jjp(fn!+i - 2/m + /m-i) + 0(A 3 ). 

3.1.2. F n —f ni ^ n2 —f(xi, x 2 ). p=2 is the lowest nontrivial value of p for which we can consider 
F n as a function of the two scales, n\ and n 2 . Taking 1=1, from eq. ( |3.7[ ) we have 

(3.23) 



p s , , df(xi,x 2 ) h 2 e 2 d 2 f(xi,x 2 ) 2 df(xi,x 2 ) 
F n ±i = f(Xi,x 2 ) ± fte 1 — ± he 



dxi 



dx\ 



dx 2 



h 2 e 



3 d df{x 1 ,x 2 ) 



+ 0(e 4 ). 



dx\ dx 2 

If F n is a slowly varying function of order two in rii, it might be of order one in n 2 . In this case, 
cq. (3.23) becomes 

F n ±i = f(xi,x 2 ) ± he 



+ 0(e 3 ). (3.24) 



df(xi,x 2 ) h 2 e 2 d 2 f(xi,x 2 ) ± ^ df(x 1 ,x 2 ) 
dx\ 2 ctef dx 2 

Moreover, from eq. (3.15) it follows that d X2 =A ri2 , <3| =(A ni ) 2 and d Xl d X2 =A ni A ri2 . Then eqs. 
( 3.23|3.24 |, written in terms of differences instead of derivatives, are given by 



Fi±l = fnx,n 2 i ^^(fn 1 + l,n 2 — fn 1 —l,n 2 ) 

2N 2 ^'^™ 1 + 1 '™ 2 _ ^^'™ii n 2 ~^ /711 — l,7i 2 ) i 2^y2 (fn 1 ,n 2 + l — /ni,n 2 -l) 
+ ^^(/iu + l,n2 + l — /m-l,n 2 + l — /rai-fl.na-l + /m-l,« 2 -l) + 0(JV 4 ) 



(3.25) 



and 



F n ±l — fn ly n 2 i \fm+l,n 2 fni — l,n 2 ) (3.26) 
"I" 2_/y2 (/»il + l:"2 — 2/ ni! ri 2 "I" /ni-l,n 2 ) i 2jy2 (^™i> n 2 + l — fn 1 ,n 2 — l) + 0(N ), 

respectively. 

3.1.3. F n ^ m —f nit . mum2 —f(xi,ti,t 2 ). In this case we have 

_df(x 1 ,t 1 ,t 2 ) 



F n ,m±l = f(xi,ti,t 2 ) ±T£- 



r 2 e 2 d 2 f{x ll t ll t 2 ) 2 df{x 1 ,t 1 ,t 2 ) 3 
± re — h C(e J 



fit? 



di 2 



and 



x , df(xi,tx,t 2 ) h 2 e 2 d 2 f(x 1 ,t 1 M) ,* 



9xi 2 

In terms of differences, the last two equations are given by 

1 
2~iV 



-^n,m±l — /ni, mi, m2 n ii r (/ni.mi+l,m2 /ni,fni- 1,7712) 



~t" 2jy 2 (f n l- m l + l-, m 2 2f nimim2 + / ni ,7711 — 1,7712 ) 

i 2jY"2 7711, 7712+1 _ All, 7711, 77l 2 — l) ^ 0(N ) 



(3.27) 



(3.28) 



(3.29) 
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and 



^n±l,m /ni, mi, 7712 (/t7-i+1,777i ,7712 /711 — 1,7711,7712) 

~t~ 27\T2 (/711 + 1, 7711, m 2 2fn 1 ,m 1 ,m 2 fn\ — 1, 7711 ,7712 ) "I" 0(N ). 



(3.30) 



For future use we can further rescale the lattice with some extra parameter by defining ni = ^jp, 
mi = and m2 = where the order 1 parameters L\ and L 2 are divisors of N and TV 2 
respectively if we require that ri\ and n 2 be integer numbers. In this case, the equations (3.271 

df(x u t u t 2 ) , T 2 L 2 e 2 d 2 f(x 1 ,t 1 ,t 2 ) 



and (3.28) become 

F n , m ±i = f(xi,h,t 2 ) ± rL 2 e 
2 df(x 1 ,t 1 ,t 2 ) 



±T€ Z 



dU 



and 



F n ±i,m = f(xi,h,h) ± hLie 



df{x u t u t 2 ) , h 2 L 2 e 2 d 2 f(x 1 ,t 1 ,t 2 



dxi 



dx\ 



Moreover, from equation (|3.31|) we have 

F n , 771-1 = T z L 2 e 



^71,771 + 1 ^F nrn 



2 , 2 2 9 2 f{xi,ti,t 2 ) 



0(e 3 



In terms of symmetric difference operators these equations can be written as 

L 2 

-^71, 771+1 /ni ,7711 ,7772 ^(/ni.TTH+l^ ./n 1 , 771 1 — 1 ,7772 ) 

^2 

"I" 2^y2 (/™l, m l+l, m 2 — ^fn 1 ,m 1 ,m 2 fn 1 ,m 1 —l,m 2 ) 



± 



2N 2 



(/ni ,77li , 7712 + 1 /ni, 7711,777.2 — l) "I" O (N ), 



(3.31) 



+ 0(e 3 ). (3.32) 



(3.33) 



(3.34) 



-^71+1,771 — /7I1 ,7711 ,7772 2 1 + ^ ' m 1 < m2 ./ni — 1,7771 ,7712 ) (3.35) 

£2 

1+1, mi ,7712 ,mi,m2 fn\ — l,7n\,rri2 

) + 0(N- 3 ) 



and 



£2 

-^71,771+1 ^F n m ~t- F n rn —l — (^7li ,7171 + 1,7712 2/77,1,777,1,7772 ~1" /ni ,77li — 1 ,7712 ) (3.36) 

+ 0(A^ 3 ). 

The last three equations will be used in the following Section to apply the multiscale method to 
the biatomic lattice model we introduced in Section 2. 

4. MULTISCALE REDUCTION OF THE DISCRETE BIATOMIC SYSTEM. 

4.1. Equations of motion. In the equations of motion of the biatomic chain (see eqs. 
( 2.2|2.3 )), the nonlinear terms (proportional to fi\ and (3 2 ) are of order e respect to the remaining 



terms, and thus we can use perturbative methods to look for approximate solutions of x n (t) and 
y n (t). This has been done in 1993 by Campa et al. [9] using the multiscale perturbative method 
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with just the lowest order differential terms. In this way, performing at the same time a multiscale 



expansion and a continuum limit they were able to reduce the system to the NLSE (5.1 ). 

Here we discretize time and look for completely discrete equations, i.e. passing from the 



differential terms in the expansion (see eqs. (3.21 3.23 3.24 3.27 3.28 3.31 3.33 1) to difference 



terms corresponding to the lowest order of slow varyness p, i.e. to eqs. (3.22 3.25 3.26 3.29 



3.30 3.34 ■ 3.36 ). To discretize time we replace the time t with a discrete variable m, so that t=mT, 



where r is the temporal scale. Thus, when r reduces to an infinitesimal quantity and m approaches 
infinity in such a way that t remains finite we recover the continuous case. We consider the simplest 
approximation of the second derivative by differences using a central difference so as to get a real 
dispersive relation. The discretized equations of motion are given by 



ITll (■^n,m+l 2x n ?n -\- X n , 
m 2 (yn,m+l - ^y n ,m + Vn,: 



-l) — ki(y nm — x nm ) — k 2 (x nm — y n -i,m) (4-1) 

~t~ ^[/^l(yn,m *£n,m) fi 2 (x n ^ m yn — l,m) ]> 
-l) = -ki(y n ,m - Xn,m) + ^(in+Lm ~ Vn,m) (4.2) 

where x Ilim =x„(mr), y n , m =y n (mT) and mi, 2 =^^- . We are looking for x n . m and y n , m as bounded 
solutions written as a modulation of the harmonic wave solutions of the linearized equations which 
one obtains when setting e — 0. The harmonic waves are given by 



E n ,r. 



i[kn— ui(k)m] 



(4.3) 



with Lu(k) real for any real value of k. The physical reason for choosing harmonic waves is that 
the atoms of the chain make only small oscillations around their equilibrium position. When we 



introduce this ansatz into equations (4.1) and (4.2), we realize at once that the solution of the 



nonlinear equations of motion can be represented as a modulated linear combination of harmonic 
functions. 



A solution of the linear part of eqs. (4.1) and (4.2) (/3 1 =/3 2 =0), written in terms of the harmonic 
waves (4.3), is given by 



where 



B _kt + k 2 + 2mi(cosw(fc) - 1) _ k 1 + k 2 e ik 

A~ ~ k% + k 2 e-' lk ~~ fci + k 2 + 2m 2 (cosw(fc) - 1) ' 



with the dispersion relation 

Ul(k) — arccos{l-j^ 



[(fei+fe2)(mi+«i2)±\/ (fei+fe2) 2 (mi+ni2) 2 — 16k 1 k 2 m 1 m2Sin 2 |] }. 



(4.4) 



(4.5) 



It can be proved that the term inside the square root of the dispersion relation is always positive, 
so that the argument of "arccos" is always real. 



In equation (4.5), the positive sign corresponds to the optical branch w opt (fc), whereas the 
negative one to the acoustical branch uj ac (k). It can be proved that the function w(fc) is real for 
all real values of k iff the temporal scale r satisfies the following inequalities: 



AM 1 M 2 



(fci +k 2 )(M 1 +M 2 



for the optical branch, and 



r< 4 



8MiM 2 



(fci + fe 2 )(Mi + M 2 ) - y/(k! + fc 2 ) 2 (Af! + M 2 ) 2 - \^k x k 2 M x M 2 



(4.6) 



(4.7) 
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E 
o 



0.5 F 



0.4 



0.3 



0.2 



0.1 



0.0 



Optica 



Acoustica 



Fig. 4.1. Graph ofoj(k) against k, with k lying in the interval [0,n]. We have chosen M\=l, Mi=1.5, 
k x =l, k 2 =0.3 and t=1(T 1/2 . 




Fig. 4.2. Graph of ui{k) against k, with k lying in the interval [0,nj. The parameters Mi, M2, k±, and 
ki are the same as in Fig. 4, but t=t . 



for the acoustical one. It is easy to show that r a is always larger than t d . In Figs. (4.1 4.2) we 
show how cj(fc) varies as a function of r. We have chosen, following Campa [9] , the following 
numerical values for the parameters, M\=l, M2—1-5, ki=l and fc2=0.3, so that r ~l. 358732 and 
r a ~2. 910816. So the obtained threshold values r Q and r a are consistent with the request that r, 
the discretization parameter, be smaller than one. 



Let us seek a finite amplitude solution of the nonlinear system (4.1 



4.2). To do so, we write 
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x„, m and y n ,m in terms of the harmonics of the linearized equation (4.3) 

oo oo 

s=0 s = l 

oo oo 



(4.8) 
(4.9) 



s=0 



where, as the variables a;„, m and 2/ n , OT are real, (G s n m , m ) are the complex conjugates of the 
modulation coefficients (G^ m , We choose G* im =ff* 1>mi)ma and H^ m =h s ni mi m2 as slowly 

varying functions of the second order in ni and mi and of the first order in m2, defined in such 
a way to avoid secular terms. Moreover we expand the functions gf ll m m2 and ft* m m2 in the 
small parameter e. So we have: 



oo 



(4.10) 
(4.11) 



z=o 



4.2. Derivation of the equations of motion. Substituting ansatz (4.8 4.9 1 into the equa- 
tions of motion (4.1 4.2) and taking into account eqs. (4.10 4.11) we get two equations of the 
form 



oo oo 



oo oo 



rr e lpW) (E Y +W e l F (s ' l) (E V 



0, 



(4.12) 



s=0 Z=0 



s=0 i=0 



where the Fn^mi ,W2 are function only of the slow variables. As (E n ^ m ) s and (E ntm ) s are indepen- 
dent functions, its coefficients must be equal to zero. So for each power of {E n>m ) and e we get sets 
of equations Fn l '% 1 ^ ri2 — for the slow varying modulation coefficients s4i',mi,m 2 and ftni,'mi,m 2 
together with their complex conjugate. 

4.2.1. e°. We look here for the linearized terms. In this case, the coefficient of the zeroth 
harmonic satisfies the equation 

(o,o) = ft (o,o) u 13) 

whereas the coefficients of the first harmonics gives a set of two equations that are identically 
satisfied when uj(k) satisfies the dispersion relation (4.5) and 

1,(1,0) 

'mi ,mi ,7712 



„(i>o) 

It can be proven easily that, for q>2, gl'fin ,m 2 = '4 ? i,m I ,m 2 =0. 
4.2.2. e 1 . The coefficients of the zeroth harmonic are 



(4.14) 



^■ 1 )(a;i > ti > t 2 ) = ff (0 ' 1 Ha:i,ti,«2) + 



fei 



(4.15) 



fcl + k 2 



\J3 2 \1 - e~ lk r\ 2 - /3i|l - r\ 2 ]\g {lfl \x ll t ll t 2 )^ 



or 



/j(0,l) 



= „(o,i) 



^lfe , (0,0) (0,0) 



2(fci + fc 2 ) 



l,mi ,m2 ' 



(4.16) 



fci 
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depending if we use the expansions in terms of derivatives or differences. 

For s=l we find a system of two equations in the two unknowns giV,^ "' "' ' ' ' 
system is compatible only if 



iti2 and /ini,mi,m2* This 



(1,0) = (1,0) 
imi ,mi ,m2 — "'"'2 ,m-2 ' 



(4.17) 



where n 2 — n\ — mi and 



fc(l>i) = 7- (i.i) 

t 'ni,mi,m2 yrii,mi,m2 



~—ik 



2z summit fc + fc 2 re~ - ( 1)0 ) _ (i,o) , 

2(fci + fc 2 e~ ife ) 1 wn J +l,mj 5n 2 -i,m 2 J> l 4 - 18 J 



where 0;^=^ =^, with L\ andL 2 given in Appendix A. 1 byeqs. (A. 2) and (A. 3). The differential 
version of eq. (4.17) is 

g^°Hx 1) t 1) t 2 )=g^(x 2 ,t 2 ), 
where x 2 =hn 2 =h(ni — m\)=x\ — ^t±, and 

,(11), . . x (ii), . ,x , 2isinwTOitJ ifc + k 2 re- lk dg(x 2 ,t 2 ) {1 - 0) , . 

^'^(xijiijia) = rg^^ 1 ^ xi,ti,i 2 + ; — -7 n: ^1 ~ ■ (4.19) 

k\ + k 2 e~' h ox 2 



For the second harmonic we get 



0(2,1) 

i/ri-i ,mi f ?7l2 

ni ,mi t m2 



iW X ' 0)2 

i,Vni jWl ,7712 1 

ifni ,mi ,77i2 ' 



(4.20) 
(4.21) 



where Ki and if 2 are given in Appendix A.l by eqs. (A. 5 A. 6). It can be easily proven that, for 

,(9,1) 



7772 "-77 1, 7771, 7772" 



=0. 



4.2.3. e 2 . Taking into account eq. (4.17), the zeroth harmonic gives a system of two equations 
that is satisfied only if 



r2/„(0,0) , (0,0) _o„(o,0) x_r c°n„(i,o) |2_|„(1,0) i 2 \ 
^ivyn 2 + l,m 2 ^ 9n 2 -l,m 2 z i/"2,m2j ~ M j Ui/n 2 + l,m 2 l li/ra 2 -l,m 2 l J 



j_ r c 1 r (1,0) /— (1,0) -{i.u) \ _^±,u; / V X , U J \\ 

+^lTU?™2> m 2U7ri 2 +l,m2 _ 9n 2 -X,m 2 ) + 3"2,m2 l5n 2 + l,m 2 — 9n 2 -l,m 2 )l' 



=(1.0) 



(1,0) / (1,0) 



,(1,0) 



(4.22) 



where c° and c 1 are two real constants given in Appendix A. 3 Defining 



-Li(a m) - o (0 '°M - — rio (1,0) l 2 + lo (1 ' 0) l 2N ) 

1 n 2 ,m 2 — -^Hi/n 2 + l,m 2 !Jn 2 ,m 2 ) q ^"a+l."^ ' ^iyn 2 ,m 2 l J 



L( (i.°) j^ 1 ' ) +o( 1 >°) o (1 < 0) ^ 

2 \3n 2 ,m, 2 h)n 2 + l,m 2 1 im 2 , 7712^712 + 1,^2 /' 



eq. (4.22) reads: 



(4.23) 



.4 



n 2 + l,m 2 



Thus A n2ir) 
a slowly varying function in n 2 we have 



(4.24) 



C(m 2 ), where C(m 2 ) is an arbitrary function of m 2 . Using the fact that gn^m2 ls 



1.(1,0) , 2 N 
\gn 2 ,m 2 \ ) 



ij Hffn 2 + l.m 2 — 9n 2 -l,m 2 ) ~ C \\9n 2 + l,m 2 

+C (ffn 2 ',m 2 3„ 2 ' +lim2 + gn 2 \m 2 g n2+ -i <m2 ) + C(m 2 ) 



Eq. (4.25) written in terms of the derivatives reads: 



(4.25) 



(4.26) 
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If we transform the derivatives of eq. (4.261 into differences (using again eq. (3.15), and recalling 
that x 2 =hn 2 ), we have 



n l\9n 2 + l,m 2 



?£°LJ = 2(c° + c 1 )^',! | 2 + C(m 2 ), 



(4.27) 



an equation simpler than eq. (4.25). This difference is due to the fact that eq. (4.26) is obtained 



using the Leibniz's rule and an integration, while in the case of eq. (4.25) the Leibniz's rule is not 
applicable as we deal with differences. 

(1 2) (1 2) 

Finally, for s=l, we get a system of two equations in the two unknowns, g)i 2 \m 2 and hn 2 \m 21 
which is compatible and not-secular only if 



zi >Hf, l2 ,m 2 + l 



(1,0) 

9n 2 ,7712 



l) + B 2-£'l(fl , „ 1 '+l, mj 

Qn 2 , 777 2 



2o (1 ' 0) \ 

£yn 2 ,m 2 ) 



+Bz\gn 2 ,m 2 
R ( (1,0) -(1,0) (1,0) (1,0) 

&5\9n 2 , 777 2 ffn2 + l, 7712 + 5772 , 7772 3n2 + 1,777 2 



-a {1 ' 0) 

2 , l-(l.O) , 2] 

77l 2 I + Iy™2,"l2 I ) 

B 6 C(m 2 )}gn 2 % 2 = 



(4.28) 



Here the coefficients Bi (i=l,...,6) are real and given in Appendix A. 3 This is a NLSE on the 



lattice. At difference from the standard discrete-time NLS equation presented by Ablowitz and 
Ladik this is completely local but not integrable 28 S! 39 . In the development of x rhm and y. 

gn 2 °m 2 is the main term which multiplies e° and E n ^ : 



77,777 c* iAV -A 1/77,777, 

If we require that gn 2 % 2 and hn 2 %i 2 are 



localized with respect to n 2 , we have to set C{m 2 ) — and eq. (4.28) becomes 



lt >l (.g„ 2 ,m 2 +i 



,(1,0) 



,(1,0) 



..(1.0) 

i>n 2 — 1,7712 



9n 2 ,m 2 -l) + ^2^1(5712 + 1,7712 



2o (1 ' 0) ) 

^yn 2 ,m 2 ) 



\1„W> >{j3(\ n W) 12 , |„U,o; |2\ 

I 5772,7772 + |-t ? 4U5„2 + l ! 777 2 I + |5t72,777 2 | j 



, R / (1,0) -(1,0) , -(1,0) (1,0) u 

+ £"5(5712, 777257l 2 + l,m 2 ' ^772 ,777 2 5t7 2 + 1 , 777 2 > J 9 



(1,0) 



= 0. 



(4.29) 



5. CONTINUUM LIMIT OF THE DISCRETE NLS. Eq. (4.29) is obtained from eqs 



(4.1 4.2 ) by discretizing the continuous time variable. This discretization was necessary to be able 
to solve the I = 1, s = 1 system which otherwise would have been an unsolvable linear differential 
difference wave equation. By discretizing we get a discrete wave equation whose general solution 
is given by an arbitrary function of a discrete variable. 



It is interesting to perform the limit when the discrete time mi is transformed into a continuous 
^-variable. To do so, we take the limit when r goes to zero and m tends to oo in such a way that 



the product rm = t is finite. So eq. (4.29) becomes the integrable NLSE 



. Ai dg^(z 2 ,t 2 ) + A fg^\z 2 M) + [Mg im {z2 ^ + Ai c(h)]g (1 >°\^t 2 ) = 0, (5.1) 

Ot 2 OZn 



where ^2=li m r-j-o li m 777^oo rm 2 arid z 2 = j,{ni~ 4rii) is a new continuous variable. The coefficients 
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Ai (i=l,...,4) in this limit are finite and real, and are given by 

, (Mi + M a )(Ai + fc 2 ) - 2M 1 M 2 fl 2 



2 



A 1 = lim 2tB 1 = -fl 

r^o k 1 + k 2 - M 2 Q 

A = lim B = t(Ml+M 2 )('=l+'°2)-A^lM2" 2 )](",fc) 2 -MlM2(n !fc ) 2 -fc 1 fc 2 cosfc 
2 r _;.o 2 fc 1 +fc 2 -M2" 2 ' 

A 3 = Iim(fl 3 + 2£ 4 + 2B 5 ) = lim(B 3 + 2(c + Cl )B 6 ) 

= -2/3 2 (J?-l){(J?-l)|J?-l| 2 2fc2(1 - c ° 3fc) p (Ml + M2) "% 2 l^|l-J?.| 2 } 



2/3|(l-fl e ' fc ){-(l-J?. e -' fc )|l-J7. e -- fc | 2 2fc i (1 - c ° afc) - (M i +M2) " + 2( ^"' -^ ll-fle"'"! 2 } 



£"2 2 _,_ 2(1 



fei+fe 2 

^|l-i?| 2 } + 2 5 A 4 , 



2(1 --Re" ;h d , 2i 



where 



fci + fc 2 

. ,. o Mall -Jte-**| 2 + Mill --RI 2 

A 4 = hmB 6 = — — . 

T-s-o fci + k 2 



2p 2 k 1 \l-Re-' k \ 2 + 2^ 1 k 2 \l-R\' 

g = lim(ci +c 2 ) 



and 



>o v 2; (Mi + M a )(fci + /c 2 )(^,fc) 2 - fcifc 2 ' 
L> = [fc 4 + fc 2 - Mifi 2 ][A:i + fc 2 - M 2 ft 2 ] - (fc 2 + k\ + 2/cifc 2 cos2fc), (5.2) 



_ .. fci + fe 2 -Mifi 2 

ft = lim r 



■ >"6 ' fci + fc 2 e~ lfe 
f2(fe)=lim T _ ) .o ^r^- gives back the continuous dispersion relation [9]. 

6. CONCLUSIONS. In this work, introducing the concepts necessary for applying the 
perturbative multiscale method to discrete equations we have obtained a rescaled discrete equation. 
We have applied this technique to a biatomic chain model. In this way we have shown that we 
can perform in a coherent way a multiscale expansion on the lattice. If we want to remain on the 
lattice and want to avoid nonlocality then we need to restrict ourselves to slow- varying functions. 
This restriction on the class of function implies that some of the properties of the starting system 
will be lost. Among them by sure that of the integrability, which is strictly related to the analytic 
properties of the solutions. 

We have found that g 1 - 1 ' ' (the slowly varying coefficient of the first harmonic) satisfies a totally 
discrete local version of the discrete NLSE. One interesting feature of our discrete NLSE is that, 
when we perform the continuous limit in the time variable, the spatial variable becomes continuous, 
and we get the continuous integrable NLSE (5.1 1 as in the work by Campa et al. (9j. 



Appendix A. Explicit Formulas. 

A.l. gi 2 'm 2 and hnimz- Let us consider the expansion of the equations of motion with 

(1,1) 

ni ,mi .m 2 , 



l=s=l. In this case we get a system of two equations in two unknowns, gn 1 'sn 1 ,m 2 and hP~'^ 



that is compatible only if 

,(1.0) 



[(fci +fc 2 )(mi +m 2 ) +4mim 2 (cosw- 1)] sin(w)L 2 (g^^ i+1 . m2 - 
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It is convenient to choose 



L\ = S , sin(cj)[(/ci + k-i){rn\ + 7712) + 4toi?ti2(cosu; — 1)] 



(A.2) 



and 



L2 = Sk\k% sin k, 



(A.3) 



where S is a real number such that L\ (X 2 ) is an integer number. In terms of L\ and L 2 
the dispersion relation becomes u)^ = j 2 -. With this choice of L\ and L2, and assuming that 

(1,0) 
9ni,roi,fnj 



ln 2 °m 2 , with ?i2=rti — mi, we find that eq. (A.l) is satisfied. Thus the system of 



equations we are studying is compatible, and leads us to the eq. (4.181 



A.2. The discrete NLSE. In this Appendix, we show the steps necessary to find the discrete 



NLSE (4.29). First, we take the equations of motion, and select the harmonic s=l with 1=2. In this 

. (12) (12).. 

way we get a system of two equations in the two unknowns gn 2 ,m 2 an d hn 2 \m 2 i which is compatible 
only if the nonhomogeneous first order difference equation 



[{ki + k 2 )(m 1 + m 2 ) + 4mim 2 (coscj - 1)] sin(w)L 2 (^ 1 i ' 1 ^ i+ljm2 



*7u,mi-l,m 2 ) + k\k 2 sm ^l(flni+l,TOi.mj 

_ t? ( (0,0) (0,0) (1,0) \ 

r I9n 2 + l,m 2 ' 9n 2 -l,m 2 i 9n 2 ,m 2 ), 



o (1,1) ) 

i?ni -l,mi ,rn,2 / 



(A.4) 



is satisfied. Here F=F(g^J 1 m2 ,gh 2 ,m 2 ) is a given nonhomogeneous term. As the l.h.s. of this 
equation is the same as that of eq. (A.l) (but with g^±im 2 re pl ace d by g^n 2 ±im 2 )i the terms 



depending on gt 1 ' ' contained in F lea d to secular terms for the unknown g^ 1 ' 1 ^ . To avoid secular 



terms, we must set F=0 and eq. (A.4) g ives g k,' 



,(l,i) 



(1,1) 

9n 2 ,m 2 ■ 



If we substitute <?(°'°) given by eq. (4.25) into F=0, then this condition will give eq. (4.29) 
written in terms of g^'°\ 



A.3. Constants. We give here the expressions of the coefficients appearing in eqs. (4.20 



4.21 4.25 4.29) 



1. Eqs. (4.20 4.21) 



Kn 



{A(r 
,02(1 -re 
{Mr 
182(1- 



- fe ) 2 [fc 1 
2 



k 2 e lk - r(h + k 2 e- 2lk )] 



k,e lk 



-ik\2\ 



k 2 e 2ik 



k. 



k,e 2lk 



-r(fc 1 e 24fe + fc 2 )]}/{r J D}, 
r(h + k 2 e- lk )] 

r(k ie 2ik + k 2 e* k )]}/{D}, 



(A.5) 
(A.6) 



where 



D = [2mi(cos2w - 1) + k\ + fc 2 ][2™ 2 (cos 2u> - 1) + ki + k 2 ] 



k\ + 2kik 2 cos2fc). 



(A.7) 



2. Eq. (4.251 



-2k, 



ik\2 



+ /3i|l 



(mi + m 2 )(ki + fc 2 )(w,fe) 2 - kik 2 ' 



2/3 2 (h 



k 2 )\l 



-ife|2 



(mi + TO 2 )(fci + k 2 ){uj.k) 2 - k\k 2 
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3. Eq. (4.29) 



B, 



sin (w) 



(mi + m 2 )(fci + fc 2 ) + 4mim 2 (cosw — 1) 



2m 2 (cosu; — 1) + k\ + fc 2 



[( mi + m2 )(fc 1 +fc 2 ) + 4 mi m 2 (cos 



-l)]co.(o;)(^ fc )- i 
k a +»2+Sm}(coeu- 



B 3 = - 2 ^(f-l){(r-l)|r-l| 2 2 ^ (1 - coafc) + 2( gi +m ^ (c ° s "- 11 

|l-r| 2 } 



2(^-1)m _, ai 



fci + fc 2 

+ 2fe 2 (l-^f e "°){-(l-re-"')|l-r e -' fc | 2 2fcl(1 - ; ° sfc)+2( g 1+ " ,2)(coa "- 1) 

2( re - ifc - 1) 



fci + fc 2 



|l-re-**| 2 } 



+ 2ff 1 fe(f-l){(f-l)(l-re-' fc ) 2 - 2e ' "■il'°-»-')-""3(""-') 

l-re~ lfc | 2 } 

-2 e - 2it m 1 (c» SM -l)-2m 2 (c.s M -l) 
' " J D 

|l-r| 2 }, 



fci + fc 2 

+ 2£ 1 /3 2 (l-fe" ! ){(r-l) 2 (l-r e lfc )- 

2(l-re- ife ). 



+ 



fci + fc 2 



£ 4 = 
B 5 = 

Ba = 



c°B 6 , 

fei/3 2 |l-re- tfc | 2 + fc 2 /3i|l-r| 2 
fci + fc 2 
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